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I. INTRODUCTION 

Molecular systems with close-lying electronic states 
possess an electronic structure that is dominated by 
strong static electron correlation. Examples are (i) 
molecules far from their equilibrium structure and (ii) 
many transition metal complexes, where static correla¬ 
tion is often sizable. A method tailored to recover static 
correlation is the complete active space (CAS) ansat JD 
for the wave function. This ansatz defines an active 
space of n act electrons in lV act orbitals, in which the 
exact solution of the electronic Schrodinger equation is 
obtained by considering a full configuration interaction 
(FCI) expansion of the wave function. Accordingly, the 
CAS wave function is defined by the number of active 
electrons and orbitals and the active space is denoted as 
CAS(n act ,A r act)- 

In the conventional CAS ansatz, the construction of a 
FCI expansion leads to a factorial scaling with respect to 
increasing values of n ac t and iV act . As a consequence, 
active orbital spaces limited to about CAS(18,18) are 
computationally feasible. By contr ast, the density ma¬ 
trix renormalization group (DMRGj-^l algorithm, origi¬ 
nally developed for interacting spin chains in solid-state 
physics, iteratively converges to the exact solution in a 
given active orbital space with polynomial rather than 
factorial coslP. With DMRG, active orbital spaces are 
accessible that are about five to six times larger than 
those of standard CASSCF algorithms. To exploit this 
advantage, an increasing number of quantum-chemistry 
DMRG implementations have emerged since the late 
1990P23. Methods, with (DMRG-SCF) and without 
(DMRG-CI) a simultaneous optimization of the orbital 
basis, were devised and a few comprehensive reviews are 
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also available® 36 . Yet, even with larger active orbital 
spaces at hand, essential parts of the remaining dynami¬ 
cal electron correlation cannot be efficiently accounted for 
within a DMRG framework. Multireference Cl (MRCI) 
in its internally contracted form 36 and mul tireference 
perturbation theory (MRPT) approaches 3 '' 3 ® have been 
combined with DMRG. Such approaches may be sum¬ 
marized as ’diagonalize-then-pertub’ in the state-specific 
case and as ’diagonalize-then-perturb-then-diagonalize’ 
in the state-average quasi-degenerate cas^®. The first 
step aims at the inclusion of static correlation effects in 
the zeroth-order Hamiltonian while capturing dynamical 
correlation in subsequent steps. Although potentially ac¬ 
curate, such a strategy comes with a caveat. This caveat 
is rooted in the need for higher-order reduced density ma¬ 
trices (n-RDMs, with n > 2) of the active space DMRG 
wave function which leads to a considerably steeper scal¬ 
ing of the MRPT or MRCI step compared to the pre¬ 
ceding DMRG(-SCF) step. A notable exception is the 
very recent development of a variational perturbation ap¬ 
proach that exploits the matrix-product structure of the 
DMRG wave function and optimizes the first-order cor¬ 
rection to the wave function iteratively by the DMRG 
protocol®. 

In this paper, we pursue a yet unexplored strategy to 
effectively treat dynamical electron correlation within the 
DMRG approach. We propose a simultaneous treatment 
of dynamic and static correlation using density functional 
theory (DFT) for the dynamical-correlation part, rather 
than aiming for a conventional two-step approach. As 
a consequence, the overall scaling cost does not exceed 
that of the DMRG calculation. This feature will be¬ 
come particularly advantageous when the approach is 
applied to large systems. In turn, the coupling of DFT 
to a CAS-type wave function method provides the re¬ 
quired flexibility for cases where static correlation be¬ 
comes important and where DFT is likely to faiP^H®. Al¬ 
though hybrid approaches between wave function theory 
(WFT) and DFT have not been considered for DMRG 
yet, they have already been studied with other wave func- 
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tion methodd^^. While approximate DFT functionals 
will introduce errors on the absolute scale, relative ener¬ 
gies for hybrid WFT DFT approaches can be obtained 
with good accuracy®^. To avoid the double-counting 
problem of electron-correlation effects, we t ake ad vantage 
of the method originally proposed by SavirP^ 47 which is 
based on a range separation of the two-electron repul¬ 
sion operator into a short-range and a long-range part. 
The short-range part of the electron interaction is then 
treated by DFT, while the long-range part is assigned 
to the WFT approach. The resulting hybrid WFT- DFT 
method is denoted ( long-range ) WFT -short-range DFT 
(in short, WFT-srDFT(P^l. In this work, we present the 
(long-range) DMRG-short-range DFT variant, abbrevi¬ 
ated as DMRG-srDFT. The theory is formulated gener¬ 
ally and applies both to traditional CAS configurational 
interaction (CAS-CI) and to DMRG. It further offers a 
natural extension to excited states. 

Besides the dynamical correlation problem, a criti¬ 
cal issue for all CAS-based methods is the composition 
of the active orbital spac^ S7 i It is a well-established 
procedur^SHSSl ma k e the orbital choice based on 
natural orbitals (NOs) and their occupation numbers 
(NOONs). In addition, significant deviations (> ±0.02) 
of NOONs from the Hartree-Fock limit of 2 (occupied 
orbitals) and 0 (virtual orbitals) have been considered as 
multireference indicators. In the framework of DMRG, 
two other measures, namel y the single-orbital entropy 12 
and the mutual information 14 61 ( have become popular f>2 . 
These entropy measures are calculated from the one- 
orbital RDM - in the latter case, also from the two- 
orbital RDM — and can be exploited to examine the 
multireference character of a wave function and/or trace 
the amount of static and dynamic electron correlation 
in an electronic wave functional Being orbital-based 
measures, these entropy measures share with NOs the 
appealing feature that they can simultaneously serve as 
indicators for multireference character as well as selec¬ 
tion criteria for the active space. With our new DMRG- 
srDFT approach, we will explore the effect of effectively 
treating dynamical correlation through the short-range 
DFT functional on entropy measures. 

The paper is organized as follows: In Section |TTJ 
we briefly outline the theory of range-separated hybrid 
methods and discuss the key steps for a combination 
of short-range DFT with long-range DMRG. Section |III| 
summarizes technical details of the implementation of 
our DMRG-srDFT approach. Computational details are 
given in Section [IV] before we proceed with the first ap¬ 
plications of DMRG-srDFT in Section |Y] 

We discuss and evaluate our approach in practical cal¬ 
culations. For H 2 0 and N 2 , we investigate the effect 
of different active spaces up to FCI for the long-range 
wave function along the symmetrical bond stretching co¬ 
ordinate. As their ground-state electronic structure be¬ 
comes increasingly multiconfigurational upon bond elon¬ 
gation, these molecules can serve as prototypical exam¬ 
ples for both singlereference and multireference methods, 


We highlight differences in entropy measures calculated 
with DMRG and DMRG-srDFT wave functions to show 
that srDFT produces a stable pattern of these measures 
that is rather independent of the size of the active orbital 
space. As a final application we investigate two ligand- 
dissociation reactions of d-block metal complexes taken 
from the WCCR10 benchmark set! 64 * 65 * , for which accu¬ 
rate experimental reference data in the gas phase at zero 
Kelvin are available. Section pVT] summarizes our findings 
and outlines future developments. 


II. THEORY 


In this paper, we generally work in Hartree atomic 
units and exploit the second-quantization formalism. 
Orbital indices p, q , r, s denote spatial general orbitals, 
i,j 7 k,l inactive (doubly occupied) orbitals, and u 7 v 7 x,y 
active (partially occupied) orbitals, thus fo llowin g the 
notation by Roos, Siegbahn, and co-workeriSS ®^. The 
electronic non-relativistic Hamiltonian then reads as 

H / ) kpqEpq ± — / ( Qpqrs&pqrs ± l 7 nm (1) 

pq pqrs 

where V nn is the nuclear repulsion potential energy opera¬ 
tor and the one- and two-electron integrals over molecular 
orbitals </>i(r) are defined as 

h pq = {(j) p (r)\h\(j) q (r)} = {p\h\q) (2) 

Qpqrs = (^p(ri)0 r (r 2 )|5(l,2)|^(ri)0 s (r 2 )) 

= (pr\g{l,2)\qs). (3) 

The operators h and £(1,2) are given in first quantiza¬ 
tion: As usual, h contains the operators for the kinetic 
energy of an electron and its interaction with all nuclei in 
the system, whereas £(1,2) is the two-electron repulsion 
operator 


5 ( 1 , 2 ) 


1 

In - r 2 |' 


( 4 ) 


The E P q and e pqrs operators are defined in terms of cre¬ 
ation and annihilation operators, 


Epq — / ( a}p a ajqfj and Gp qrs — Ep q E rs Ep S S qr . (5) 

a 

Then, the electronic energy for an electronic wave func¬ 
tion T can be written in terms of the one- and two- 
electron RDMs, 


D pq = {^\E pq \^) (6) 

and 

Ppqrs — (^| &pqrs 


( 7 ) 
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respectively, as 

E = (’F| J ff|'F) = E hpqD pq + ^ E QpqrsPpqrs + Kn, 

pq pqrs 

( 8 ) 

where the operator V nn is written as a potential energy 
Km of the nuclear framework since it does not depend on 
electronic coordinates, which are the dynamical variables 
integrated out in the energy expectation value. In CAS- 
type methods, parts of the RDMs will be associated with 
inactive electrons and the computational evaluation of 
the energy expression can be further simplified by split¬ 
ting it into separate contributions for inactive and active 
electrons. The corresponding formalism will be elabo¬ 
rated in the following subsection. 


A. Complete-Active-Space Configuration Interaction 


By dividing the one-electron RDM into an inactive 
part (T), D l = {D \-} = {2 Sij}, and an active part (’A’), 
D a = {D a , }, we may write the CAS-CI energy expres¬ 
sion as a sum of an inactive energy ( E \) and an active 
energy (E A ): 


^CAS-Ci = Ei + E A , (9) 

where 

Ei = -E(^3 + fij ) E\j + V nn 

= + fii) + (10) 

i 

E A = fuyE u v + 2 9uvxyP u vxy (H) 

uv uvxy 

The inactive energy E\ is equal to the Hartree-Fock en¬ 
ergy expression for the doubly-occupied orbitals. The 
matrix element f}„, denotes an element of the inactive 

j pq 

Fock matrix 

fpq = hpq + y X ^Qpqkk — 9pkqk) i (12) 

fc 

which has been defined according to Eq. (15a) in Ref. (BT 
(see also Ref. EH for explicit derivations). In E A , the 
use of the inactive Fock matrix, f l = {/£„}, instead of 
the one-electron matrix, h = {h uv }i accounts for the 
screening of the nuclei by the inactive electrons. 

The energy expressions in Eqs. © @ hold for any 
CAS-type method, including those with orbital optimiza¬ 
tion such as the Complete-Active-Space Self-Consistent- 
Field (CASSCF) method. 


B. Range-separated CAS-CI hybrids with DFT 


The two-electron repulsion operator can be separated 
into a long-range fir’) part and a short-range (’sr’) 

paxPEEtea, 


S(l,2)=^’ lr (l,2)+^’ sr (l,2), 


(13) 


involving a range-separation parameter /i. This de¬ 
composition of the electron-electron interaction opera¬ 
tor has been applied in various WFT-srDFT hybrid 

In this 

paper, the long-range 
and short-range parts of the interaction operator are sep¬ 
arated by virtue of the error functiorP^l, 


^,h Vl r>\ erf(/x n — r 2 ) 

9 (1,2) = -;-;- 

Ti - r 2 


(14) 


y,(l, 2)= 1 - e ^ | -’- a | ) , (15) 

In the following, the long-range and short-range two- 
electron integrals, g pqrs and g p T qrs , are the integrals in 
which g( 1,2) of Eq. ([3| has been replaced by <7 M,lr (l,2) 
and g M,sr ( 1,2), respectively. All two-electron integrals 
depend on the range-separation parameter /i to be fixed 
prior to a calculation. For the sake of brevity, we refrain 
from denoting this explicit dependency for the integrals 
and for the energy in what follows. 

In the next step, the short-range part of the electron- 
electron interaction energy is described by DFT with a 
tailored (short-range) functional E|j r xc [p] of the total elec¬ 
tron density 

p(r) = (’F|/3|'I') = y^(\E'|Dp 9 A P9 |\E') = Q pq D pq , 

pq pq 

with n pq (r) = <p*(r)(/)q(r). (16) 

We note that the limits g = 0 and g -A oo then corre¬ 
spond to Kohn-Sham DFT and ab initio WFT, respec¬ 
tively. 

The srDFT functional is partitioned as usual in DFT 
methodology into a Hartree (Coulomb) term, and 

an exchange-correlation (xc) contribution, U x (.[p], 

^H r xc[p]=n r W+^xcM, (17) 


where 


Em = \ E D V1 9pqrs Drs = lj2 im D P0’ ( 18 ) 


pq,rs 


pq 


while the explicit form of E^ c [p] depends on the choice of 
the approximate functional. We have implicitly defined 
the short-range two-electron Coulomb potentials j pq in 

Eq. pb. 


It must be stressed that results for a range-separated 
hybrid WFT DFT approach in practice will be g- 
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dependent due to the approximate nature of the short- 
range functionals available. The same is true for range- 
separated Kohn-Sham DFT. Calibration studies for the 
latter suggest that values in the interval 0.33 a.u. < p, < 
0.5 a.u. are opt imaf 73 Studies using a CAS-srDFT 
hybrid® have shown that p = 0.4 a.u. is a good compro¬ 
mise which optimizes the amount of static correlation re¬ 
covered by the wave function part. An alternative which 
defined the optimal /r-value as the value that provides 
the lowest CAS-srDFT energy was also explored®, but 
was found to be system dependent (due to the approx¬ 
imate srDFT functional). It is therefore not surprising 
that for some systems the lowest energy can be obtained 
in either pure CASSCF (p = oo) or pure DFT (p = 0) 
calculations. 


approximation to the energy change, 

-££>*] (24) 

»/ Wf ^ Ap(r) dr 


where 


-ref.E 

Jpg 


=(^b'xc[p ref ] 


X Z7' sr 

*) = (<t>r\j^r ) lP let m (25) 


=E 


sr jjief 
jpqrs rs 5 


For finite p, the CAS-CI energy expression of Eq. ^ 
becomes 


-®CAS-Ci — E{ + Ea + -E'h [p\ + ^xc [p\ i 


(19) 


where the first two terms are identical to Eq. © except 
that all regular two-electron integrals have been replaced 
by the long-range two-electron integrals, that is, g pq . 


„lr 


g pqrs . Accordingly, the inactive Fock matrix in Eq. (12) 
is modified to 


are the matrix elements of the short-range Coulomb op¬ 
erator, and 

X IT'S 1 

<’ sr = (bplCb ref ]|b 9 ) = (<M^[p ref ]l<U (26) 

are the matrix elements of the short-range exchange- 
correlation potential. We can now define the state- 
averaged CAS-CI-srDFT method for M electronic states, 
using the reference density matrix elements 


fpq hpq A Y,( 2 %qkk Ppkqk) * 


( 20 ) 


D lei = — 
pq M 


M 

E W iE>pq, 


(27) 


One notes that this CAS-CI-srDFT energy expression 
is not linear in the one- and two-electron density ma¬ 
trices as in standard CAS-CI, because the Hartree and 
exchange-correlation terms in Eq. (19) are non-linear in 


the one-electron density matrix. We illustrate this by 
considering a linear deviation A D pq = D pq — D p e q from 
some (fixed) reference density matrix, D rei = {D p ° q }. 
The one-electron density matrix elements are thus 


n — n ref _i_ a D 

pq pq ' 1 -iU. 


pqi 


( 21 ) 


which by insertion in Eq. (16) leads to 


p = p ret + Ap, 


( 22 ) 


in an obvious notation. As E^ xc [p] is non-linear in the 
one-electron density matrix, we note that 

£f xc [p ref +Ap] ^ £?L c [p ref ] + Eg xc [Ap\. (23) 


This has the consequence that an exact CAS-CI-srDFT 
expression is state specific , and we cannot diagonalize a 
matrix to obtain exact CAS-CI-srDFT electronic ener¬ 
gies of several roots. As holds in general for state-specific 
methods, this implies that the Cl expansions for different 
states of the same symmetry will be non-orthogonal. 


However, following Pederserl®, we can define a lin¬ 
ear model providing orthogonal Cl states in the spirit 
of state-averaged CASSCF by using the following linear 


and thus p Tei = jj w iPii where the weights Wi add 
up to one. Although the Cl expansions now will be or¬ 
thogonal, the equations are still non-linear. The most 
straight forward optimization procedure is an iterative 
method, in many ways similar to Hartree-Fock theory; 
this will be described in Sec. Mil 

We proceed by noting that in any CAS-CI model 
A D pq = AD^ v because the inactive part, D\ :) . is fixed 
by definition. The state-averaged CAS-CI-srDFT energy 
for the i = 1, M selected roots can then be written as 

^CAS-CI -&1 + i %xc[ P \ 

+E(j- sr +<' f,sr ) AD -- ( 28 ) 

uv 

If only one root is used (M = 1), the formalism will co¬ 
incide with a state-specific optimization. In the DMRG- 
srDFT variant, all equations above hold, but a DMRG 
protocol (see the following subsection) optimizes the ac¬ 
tive density matrices. 


C. The DMRG ansatz and correlation measures 

A given electronic wave function can be expanded in 
terms of occupation number vectors, 

T= E ■ ■ a L ). (29) 

cri—<T L 
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Here, the configuration-interaction expansion coefficients 
are written as a coefficient tensor , i/j ITl ...< TL , according 
to the tensorial construction of the -dimensional 
Hilbert space from N act spatial orbitals. In DMRG ter¬ 
minology, each spatial orbital defines a site with four 
possible one-electron states, 


aj = {K’ ac )> l«)> l/?>, l«/3)}, 


(30) 


The quantum-chemical DMRG approach builds up the 
CAS wave function by first arranging the set of (active) 
orbitals {</>£} in a linear order according to some opti¬ 
mization recipe (e.g., according to the single-orbital en¬ 
tropies calculated from a few DMRG sweeps). In our 
second-generation, i.e., matrix-product -opera tor-based 
implementation of the DMRG algorithirpSGS eac h site 
has an associated set of operators in matrix representa¬ 
tion. While optimizing the site matrices iteratively, the 
DMRG protocol constitutes a variational optimization 
with respect to the total electronic energy. One- and 
two-electron RDMs as in Eqs. © and ([7]), can then be 
evaluated. For our DMRG-srDFT method, these density 
matrices are required to evaluate the energy expression 
in Eq. ([28]). 


In order to estimate the multireference character of 
the target molecule in terms of orbital-based descriptors, 
we exploit the fact that the DMRG wave function can 
be easily partitioned into two (open) quantum systems 
within the DMRG algorithm: One or two orbitals are 
embedded into all remaining orbitals of the active space. 
If | n) denotes the states defined on this single orbital (or 
on the two orbitals, respectively) and | j) those defined on 
the remaining orbitals of the CAS, the partitioning yields 
an RDM operator for the states defined on the embedded 
orbital(s), 


Pnn' =Y1 b')l n >( n 'IO' , |> ( 31 ) 

ij' 


where the environment states \j) and \j') are traced out. 
The one-orbital (or two-orbital) RDM evaluated as an 
expectation value of the RDM operator in Eq. (31), 


Pun' = (I'l/wW, 


(32) 


is then diagonalized to obtain four (or sixteen) eigenval¬ 
ues n ajU (or n a ^ uv ) for orbital (or for orbitals and 
0^). The one-orbital entropy 1 -^, 


4 

$u ^ ) ^a,u hl(?T aiW ), (33) 

a=l 

measures the degree of entanglement of the four pos¬ 
sible states defined on orbital u with all states defined 
on the environment orbitals. Similarly to Eq. ( |33| ), the 
entanglement of all 16 states defined on two orbitals em¬ 
bedded in the complementary orbital space of the CAS 


1. Start-lip run (if iter == 1) 

la. Calculate lr/sr integrals and {f I,,r } 

lb. Run DMRG - get {D A } and {P A } 

lc. Calculate E ( and the first energy, Eq. (11) 


2. Regular run (if iter >1): 



2a. Calculate E^, Eq. (37), from sr integrals, 
{D A " f } as well as 

2b. Use {f 1 ’ 1 '} and to calculate (update) y 

E= { f Ur D'} + E )* 


if (damping == .true.) 

{D *r„, j 

3. Run DMRG with modified one-electron integrals: set iter = iter + 1 

3a. Extract density matrices {D A } and {P A } 

3b. Calculate total energy; Eqs. (9), (37) and (39) 


FIG. 1. Flow chart for the DMRG-srDFT implementation re¬ 
ported in this work. Each iter denotes a macro-iteration, and 
the DMRG (or Cl) uses sufficient sweeps (or micro-iterations) 
to achieve convergence within each macro-iteration. 


is measured by the two-orbital entropy, 
16 


$uv — ^ ' rL a uv ln(^ Q;W ). 


(34) 


a = 1 


Eq. (34) also contains single-orbital contributions, which 
can be eliminated by subtracting the single-orbital en¬ 
tropies of or bitals u and u, which defines the mutual 
information 12 *^^, 


ffiv — 2 ($uv $u $v) (1 dyv ) - (35) 


III. IMPLEMENTATION 


Our DMRG-srDFT i mple mentation is based on an ex¬ 
isting CI-srDFT program' 8 and provides an interface be¬ 
tween a development version of the D ALTOisP^lSil program 
and the Maquis quantum chemical DMRG progranPD. 
Maquis is a genuine DMRG program based on ma¬ 
trix product states (MPS) and matrix product operators 
(MPO). Both regular and short-range integrals are cal¬ 
culated with Dalton. To be closer to the operational 
expressions used in the program we split the total energy 
in Eq. (28) into inactive and active parts as in Eq. ©• 
The inactive part then becomes 


^isrDFT _ 


ZT’i 


SA 
fixed ’ 


(36) 


where the fixed, state-average (’SA’) energy, -Ef^ d , con¬ 
tains short-range terms that depend only on the inactive 
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and the reference density matrices 


IV. COMPUTATIONAL DETAILS 


pSA _ ST' -I,sr n I 
fixed 2 / j •'i'j 2 


jrd ^ \ -A,ref.sr pv 

D ij ~ ^ Z^Juv D 
u 

1 rnsr r refi \ ^ „ ref.sr 7-v 

~^~^yLc\P J / j ^xc 


A, ref 
uv 


■^A,ref 
''uv 5 


(37) 


with the state-averaged active reference density matrix 

M 


7-\A,ref _ 1 \ ' 7 -)A,« 

^uv . ^i^uv ■ 


M 


(38) 


i=l 


Recall that the active reference density matrix is kept 
fixed during the Cl or DMRG optimization, and can 
therefore be assigned to the inactive energy. It changes 
only in each macro-iteration as it is obtained from the 
wave function of the previous iteration. The active part 
of the energy is 


^ DFT ’ l = ^+Ea 


I,sr _i_ -A,ref,sr 
uv ' Juv 


ref, S i\ D A,i 
^XC ,uv ) uv 


— fl.lr I „-I,sr , A,ref.sr , .ref.sr \ ry 

/ J \ J uv ' Juv ' Juv ' u xc,uvJ J - y uv 

uv 


uv 

/*I,lr 1 -I,s 


■)A ,i 


+5 E 


If pA,i 
Uuvxy uvxy 


(39) 


We can now define the iterative CAS-CI-srDFT or 
DMRG-srDFT procedure as follows: 


1. select an initial reference density; 

2. solve for the M Cl vectors / DMRG states used in 
the averaging with fixed reference density; 


3. calculate new reference density and energy; if not 
converged, go back to step 2. 

Figure [l] summarizes the workflow in steps 1-3 of our 
implementation. Note that integrals with four active in¬ 
dices are constant during the iterations. Thus, they are 
calculated (and transformed to the MO basis) only once. 
Moreover, E\ r is also a constant. The only quantities 
that need to be recalculated in each (macro)-iteration 
are j£’ ref ’ sr and <£,. 

We emphasize that steps 1 -3 and the scheme in Figure 
[T] hold for both DMRG-srDFT and general CI-srDFT 
(if DMRG is replaced by Cl). The first-order optimizer 
in the previous work 78 did not consider any convergence 
acceleration or damping schemes and convergence prob¬ 
lems were frequently observed. In order to (partially) 
solve the latter issue we introduce a simple, dynamical 
damping scheme. In iteration iter, we modify _D A,ref ac¬ 
cording to 

Df^ p = aD^f + (1 - a)D? t £L lt (40) 

where a is a dynamically adjusted damping factor. 


Calculations for H 2 0 and N 2 employed a Dunning 
cc-pVDZ basis sell® for O, H, and N. The DMRG- 
srDFT calculations were in all cases performed with 
the srPBE functional from Ref. [253 i.e. with the Heyd- 
Scuseria-ErnzerhoP^ short-range exchange functional to¬ 
gether with the rational interpolation srPBE correla¬ 
tion functional defined in Ref. (25] (SRCPBERI in DAL¬ 
TON). The range-separation parameter was always set 
to iJ = 0.4 a.u.H 

The structures considered in this work correspond to 
the ones used by Olsen et al P® for H 2 0 and by Chan et 
alW^ for N 2 . For H 2 0, we exclusively consider the sym¬ 
metric stretch coordinate. The truncated active orbital 
spaces for H 2 0 and N 2 comprise all valence electrons and 
orbitals required for a balanced description of the valence 
electronic structure along the stretching mode, that is 
DMRG(8,8) for H 2 0 (corresponding to ’CASB’ in Ref. 
152i ) and DMRG(6,6) for N 2 . H 2 0 was calculated in C 2 V 
symmetry and the CAS(8,8) space includes 4 orbitals in 
Ai symmetry, and 2 orbitals in Bi and B 2 symmetries, re¬ 
spectively. N 2 was calculated in the D 2 h subgroup. Thus, 
the used CAS(6,6) space correspond to 1 orbital in sym¬ 
metries A g , B 3 U , B 2 U , Bi u , B 2g and B 3 g . Active spaces 
corresponding to FCI (DMRG-FCI) within the cc-pVDZ 
basis set are DMRG(10,24) for H 2 0 and DMRG(14,28) 
for N 2 . The number of renormalized DMRG block states 
m was set to m = 512 for H 2 0 while for N 2 higher m val¬ 
ues of up to m = 2048 were required for technical reasons 
to achieve similar convergence as in Ref. (55] Accordingly, 
we specify the DMRG data as DMRG(n ac t,A r ac t)[TO]. All 
calculations have been performed as state-specific ones, 
i.e., with one root. 

For the ligand-dissociation reactions, we applied the 
structural models depicted in Figure [2j which were trun¬ 
cated compar ed to the original metal complexes in our 
previous wor ]j£iMI on the WCCR10 benchmark set (es¬ 
sentially, large mesityl and aromatic groups were replaced 
by methyl residues). The structures were optimized with 
the BP86 functional® and a def2-TZVI®- basis set along 
with the corresponding basis set for Coulomb fitting®. 
After the structure optimization, the Cu-N/Pt-N bonds 
are stretched to 7 A in order to mimic the final prod¬ 
ucts in a supermolecular calculation for which the active 
space can be chosen in complete analogy to the opti¬ 
mized reactant complex (cf. Figure [ 2 ]). The stretched 
structures were re-optimized with fixed Cu-N and Pt-N 
internuclear distances, respectively. All these prepara¬ 
tory calculations were carried out with the Turbomole 
program (version 6.5)—1 Then, pure DFT calculations 
were carried out with the PBE functional® and the def- 
TZVP basis selpS to understand the effect of the struc¬ 
tural truncation as well as the smaller basis set compared 
to the original work!®. 

In the DMRG and DMRG-srPBE calculations, the ac¬ 
tive spaces were chosen to encompass both d-type or¬ 
bitals and ligand orbitals important for the dissociation 
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Reaction 1 



FIG. 2. Ligand-binding energies defined by ligand dissocia¬ 
tion reactions for two model reactions out of the WCCR10 
set of ligand-binding energies. 


(see Figures in the Supporting Informatioi®). The cho¬ 
sen orbital spaces are in all cases very similar between 
DMRG and DMRG-srPBE. DMRG and DMRG-srPBE 
calculations were then performed with the def-TZVP ba¬ 
sis set using the same setup as for H 2 0 and N 2 . The 
number of renormalized states was m = 2000 in all cases. 
We also investigated the dissociation energy of reactions 
1 and 2, using a larger number of DMRG sweeps and also 
lower number of renormalized states. These tests showed 
that the setup just described was sufficient to achieve 
convergence within the given significant digits (see Ta¬ 
bles in the Supporting Information^). We observe in 
general that the number of renormalized states required 
for a converged result is smaller for DMRG-srPBE than 
for regular DMRG calculations. For the two dissociation 
reactions we have also investigated the effect of using a 
different short-range DFT functional, namely the Goll 
Werner-Stoll 71 combined correlation and exchange func¬ 
tional. This combination will be denoted srPBE(GWS). 

All DMRG and DMRG-srDFT calculations in this pa¬ 
per were carried out starting from HF and HF-srDFT or¬ 
bitals, respectively. During the warm-up sweep the state 
corresponding to the HF determinant was explicitly en¬ 
coded in the MPS. 


V. RESULTS AND DISCUSSION 

A. The effect of truncating the active orbital space 

The first part of this section is concerned with the effect 
of a truncation of the CAS both in standard DMRG and 
DMRG-srDFT calculations on H 2 0 and N 2 . Total elec¬ 
tronic energies are reported and discussed in the Support¬ 
ing Information^. While these energies show that the ef¬ 
fect of truncation of the active space is smaller in DMRG- 
srPBE than in standard DMRG calculations (due to a 
’regularizing’ effect of the srDFT part on the CAS), the 
explicit energy data is not meaningful as our current im¬ 
plementation does not support spin-unrestricted srDFT 


calculations. As a consequence, the energies of the open- 
shell reaction products will be asymptotically unreliable 
(and in fact, worse than those obtained from unrestricted 
Kohn Sham DFT calculations). Moreover, the small ac¬ 
tive spaces chosen for these two molecules are already 
so large that good agreement with the FCI reference is 
obtained in the large-CAS DMRG calculations. There¬ 
fore, we postpone a discussion of energies to the next 
subsection. Here, we continue to explore the ’regulariza¬ 
tion’ effect of srDFT on the active space of the DMRG 
calculations by investigating the entanglement measures. 

In a recent stud}®, we suggested that single-orbital 
entropies, s u , and mutual information, I uv , can serve as 
descriptors to trace and classify multireference character. 
We found that a balanced active space should comprise 
all orbitals with s u > 0.5 and/or all orbitals within a 
range of 0.01 < I uv < 0.1. As these descriptors can be 
a useful guide to construct optimal active orbital spaces 
for DMRG and standard CAS-type calculation®, they 
may complement a selection procedure based on natural 
orbital occupation numbers. The single-orbital entropies 
and mutual information for H 2 0 and N 2 calculated at 
various stretched structures are summarized in Figures 
[3] and [fj respectively. In these entanglement plots the 
magnitude of the single-orbital entropy for each orbital 
is determined by the size of the corresponding red circle 
while the mutual information is encoded by line color and 
thickness; the thicker and darker the connecting line be¬ 
tween two orbitals the larger is their mutual information. 
Considering first the standard DMRG data in the upper 
parts of Figures [3] and |4j our entanglement data confirm 
the chemically intuitive choice that all orbitals involved 
in O—H or N—N bonding need to be included in a mini¬ 
mal CAS. In addition, in accordance with the boundaries 
of the entanglement measures defined above for standard 
DMRG®, the stretched H 2 0 and N 2 molecules display 
significant multireference character. In contrast, both the 
single-orbital entropies and the mutual information are 
significantly smaller for DMRG-srDFT than for DMRG. 
For both molecules a major part of the dynamical corre¬ 
lation is indeed treated by the srDFT functional, while 
static correlation is efficiently taken care of by the long- 
range wave function. Upon bond stretching, static cor¬ 
relation becomes more dominant, but the amount of dy¬ 
namic correlation assigned to the long-range wave func¬ 
tion remains effectively constant (cf. the patterns of 
green lines for the DMRG-srPBE entries in Figures [3] 
and [4] when proceeding from left to right). 

Figures [3] and [4] further show the effect of truncating 
the active spaces. For regular DMRG the static corre¬ 
lation increases and appears to be overestimated in the 
elongated systems for the truncated active spaces. Hence, 
the effect of active-space truncation is much smaller in 
DMRG-srPBE. Therefore, srDFT has a regularizing ef¬ 
fect on the entanglement of orbitals in the active space 
as the qualitative picture provided by the entanglement 
measures does hardly change in DMRG-srDFT calcula¬ 
tions when the active space is reduced. 






3 




o 

o 







□ 

7 ) 

0 

A 

UJ 

m 


FIG. 3. Single-orbital entropies (encoded in the size of the red circles) and mutual information (encoded in the color and 
strength of the connecting lines) for H 2 0 at three bond distances. Top: Entanglement plots from DMRG(10,24)[512] with 
DMRG(8,8)[512] given below. Bottom: Entanglement plots from DMRG(10,24)[512]-srPBE with DMRG(8,8)[512]-srPBE 
given below. Active orbital labels are in black and inactive/secondary orbital labels are in gray. The numbering of orbitals is 
the same for large and small active spaces. 


From the above discussion it is also clear that the rec- measure for the CASSCF-srDFT 70 hybrid approach, 
ommended boundaries for assessing a minimum active 
orbital space with entanglement measures need to be re¬ 
vised for range-separated hybrid methods. Similar con¬ 
clusions have been drawn with respect to the boundaries 
for natural occupation numbers as active-orbital space 
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FIG. 4. Single-orbital entropies (encoded in the size of the red circles) and mutual information (encoded in the color and 
strength of the connecting lines) for N 2 at three bond distances. Top: Entanglement plots from DMRG(14,28)[2048] with 
DMRG(6,6)[2048] given below. Bottom: Entanglement plots from DMRG(14,28)[1024]-srPBE with DMRG(6,6)[1024]-srPBE 
given below. Active orbital labels are in black and inactive/secondary orbital labels are in gray. The numbering of orbitals is 
the same for large and small active spaces. 


B. Ligand binding energies in transition metal complexes 

The ligand-binding energies of the WCCR10 set of lig¬ 
and dissociation reactions are hard to reproduce by DFT 
and the PBE functional turned out to be the pure den¬ 
sity functional with the smallest overall error^J. For this 


reason, we chose two reactions from this test set to in¬ 
vestigate the potential of our DMRG-srDFT approach 
and employed the structural models depicted in Figure 
[2] Tables |T] and |TT] provide all data obtained for these 
two reactions. We will focus on the pure electronic con¬ 
tribution to the dissociation energy, i.e. to D e , as the 
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Dq data are all obtained by the same constant shift and 
are reported only as they represent the true experimental 
observables at zero Kelvin. 

First of all, we should discuss the effect that the re¬ 
duced structural model and the reduced basis-set size 
have on the dissociation energies reported in Tables [I] 
and El This assessment can be done based on the PBE 
results. As can be seen from the Tables, for reaction 1 
we increase the electronic contribution to the dissocia¬ 
tion energy D e from 247.5 kJ/mol to 257.5 kJ/mol by 
switching from the quadruple-zeta basis set of Ref. [Ml to 
triple-zeta basis set employed in this work. This value 
is then reduced to 240.2 kJ/mol by reducing the size of 
structural model. Hence, when considering a zero-point 
vibrational energy correction (see Tablejl]) the experimen¬ 
tal reference result of 218.2 is enlarged to 226.7 to obtain 
the D e reference result, which is to be corrected for the re¬ 
duced basis set by 257.5—247.5=+10.0 kJ/mol and then 
for the reduced structural model by 240.2—257.5=—17.3 
kJ/mol. The final adjusted reference value for D e is ob¬ 
tained as 226.7+10.0—17.3= 219.4 kJ/mol. In a simi¬ 
lar manner, we can adjust the D e reference energy of 
109.9 kJ/mol in Table |n| by +7.5 kJ/mol for the model- 
structure error and +12.6 for the basis-set error to finally 
yield 109.9+7.5+12.6=130.0 kJ/mol. 

In Table [T] for reaction 1, we note that the dissoci¬ 
ation energy strongly depends on the size of the ac¬ 
tive space in the pure DMRG calculations; the disso¬ 
ciation energy is increased from 132.8 kJ/mol for the 
small CAS to 173.5 kJ/mol for the largest CAS. This 
dramatic spread is not seen in the DMRG-srDFT calcu¬ 
lations, where it ranges only from 216.5 kJ/mol to 225.1 
kJ/mol. Moreover, we note that these latter energies 
are in excellent agreement with the reference energy of 
219.4 kJ/mol. Hence, the dynamic correlation captured 
in the srDFT part allows us to apply a much smaller ac¬ 
tive space and yields results in much better agreement 
with the reference energy. The agreement is also bet¬ 
ter than the one obtained within pure DFT calculations. 
Table [I] also allows for a comparison of different srDFT 
functionals. From the dissociation energy of 246.5 ob¬ 
tained with DMRG[2000] (30,22)-srPBE(GWS) and com¬ 
pared to 225.1 obtained with DMRG[2000](30,22)-srPBE 
we understand that the effect of the approximate func¬ 
tional can be larger than 20 kJ/mol and thus further 
away from the reference result (but comparable with the 
pure DFT result). 

For reaction 2, we found an even more pronounced de¬ 
pendence of the dissociation energy on the size of the 
active space for the pure DMRG data reported in Table 
|n} compare, for instance, 65.3 kJ/mol for DMRG(8,8) 
to 34.0 kJ/mol for DMRG(22,20). As this change in 
energy also increases the deviation from the experimen¬ 
tal reference energy although the CAS was enlarged and 
should have improved on the small-CAS result, it calls 
for a systematic investigation (data reported also in Ta¬ 
ble [TT| see also the total electronic energies reported in 
the Supporting InformatiorP^, which clearly show that 


the change in dissociation energy is brought about by 
the product structure only). We see that the dissociation 
energies obtained for the small- to medium-sized active 
spaces in Table |TT| only change moderately. However, for 
DMRG(18,18) we observe a significant lowering of the 
dissociation energy, which can be explained by closer in¬ 
vestigation of the orbitals that were absent in the smaller- 
CAS calculations (for orbital diagrams and natural oc¬ 
cupation numbers see the Supporting InformatiorP^) : In 
the product structure, the newly added occupied orbital 
(orbital no. 56) differs with an occupation number of 
1.966 significantly from 2.000. It is also different from 
the occupation number of 1.992 obtained for its corre¬ 
sponding orbital in the reactant structure (orbital no. 
61). This large change demands the inclusion of the or¬ 
bital pair for a correct description of the reaction, espe¬ 
cially when dynamical correlation is not accounted for 
otherwise. Notably, the DMRG(8,8) to DMRG(16,16) 
calculations do not include this orbital in the CAS of the 
product structure, and the better correspondence with 
experiment must therefore be considered fortuitous. In 
the largest CAS considered, we obtained 44.3 kJ/mol 
for the electronic contribution to the dissociation energy. 
Compared to the 34.0 kJ/mol for DMRG(22,20) we thus 
again start to improve the correspondence to the experi¬ 
mental dissociation energy, although the deviation is still 
dramatic, which might be taken as an indication that dy¬ 
namic correlation is pivotal in the dissociated fragments. 

We carried out a similar investigation of the orbitals 
for the DMRG-srPBE calculations. In this case, the 
effective inclusion of dynamical correlation renders the 
change in occupation numbers about an order of magni¬ 
tude smaller. For the orbital pair discussed above, the 
occupation numbers change from 1.996 (orbital no. 59) 
in the reactant to 1.991 (orbital no. 55) in the prod¬ 
uct. Accordingly, the dissociation energy is less affected 
(the results being 80.6 kJ/mol for the smallest CAS and 
81.7 kJ/mol for the largest CAS). We note that these re¬ 
sults are similar to the corresponding PBE result of 86.4 
kJ/mol, but still largely deviate from the reference value 
of 130.0 kJ/mol. In order to understand whether the 
approximate nature of the short-range functional can ac¬ 
count for this deviation, we investigated the alternative 
srPBE(GWS) functional which yields an dissociation en¬ 
ergy of 101.8 kJ/mol that is considerably closer the ref¬ 
erence. Clearly, changing the srDFT functional is not 
a universal solution as the srPBE(GWS) functional for 
the dissociation energy of reaction 1 yields 246.5 kJ/mol, 
thus overestimating reaction-1 reference value of 219.4 
kJ/mol. Hence, this srDFT-functional study emphasizes 
that the available short-range functionals should be im¬ 
proved, as has also been noted by otherP^. 


VI. CONCLUSIONS AND OUTLOOK 

We have presented the development and first imple¬ 
mentation of a hybrid approach that couples DMRG and 
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TABLE I. Calculated dissociation energies in kj/mol for reaction 1 obtained with the def-TZVP basis set. Do is the zero-point 
vibrational-energy corrected employing a value of 8.5 kj/mol for the zero-point vibrational energy obtained for the full complex 

with DFT(BP86)/def2-QZVPP (Refs. [MMI- 


Method 

D e (kJ/mol) 

D 0 (kJ/mol) 

DMRG [2000] (30,22) 

173.5 

165.1 

DMRG [2000] (20,18) 

169.9 

161.5 

DMRG[2000] (10,10) 

132.8 

124.3 

DMRG [2000] (30,22)-srPBE(GWS) 

246.5 

238.1 

DMRG [2000] (30,22)-srPBE 

225.1 

216.6 

DMRG [2000] (20,18)-srPBE 

227.9 

219.4 

DMRG [2000] (10,10)-srPBE 

216.5 

208.0 

PBE 

240.2 

231.8 

PBE (full complex/def2-TZVP) 

257.5 

249.0 

PBE (full complex/def2-QZVPP from Ref. IMI) 

247.5 

239.0 

Exp. (Ref. IM1) 

226.7 

218.2 

TABLE II. Dissociation energies in kj/mol for reaction 2 obtained with def-TZVP basis set. 

Do is the zero-point vibrational- 

energy corrected result employing a value of 7.3 kj/mol for the 

zero-point vibrational energy obtained for the full complex with 

DFT(BP86)/def2-QZVPP (Ref. 64 65f. 



Method 

D e (kj/mol) 

D 0 (kj/mol) 

DMRG [2000] (24,24) 

44.3 

36.9 

DMRG [2000] (22,20) 

34.0 

26.7 

DMRG [2000] (18,18) 

37.6 

30.2 

DMRG[2000] (16,16) 

56.6 

49.3 

DMRG [2000] (14,14) 

55.2 

47.8 

DMRG [2000] (10,10) 

60.4 

53.1 

DMRG [2000] (8,8) 

65.3 

58.0 

DMRG[2000](22,20)-srPBE(GWS) 

101.8 

94.5 

DMRG [2000] (22,20)-srPBE 

81.7 

74.3 

DMRG[2000](8,8)-srPBE 

80.6 

73.2 

PBE 

86.4 

79.1 

PBE (full complex / def 2-T Z VP) 

73.8 

66.5 

PBE (full complex/def2-QZVPP from Ref. IMI) 

66.3 

59.0 

Exp. (Ref. (Ml 

109.9 

102.6 


DFT using a range-separation ansatz. To analyze the 
DMRG-srDFT approach, we considered the (symmetri¬ 
cally) dissociating H 2 0 and N 2 as well as ligand-binding 
energies in transition-metal complexes. 

Although total electronic energies from small-CAS 
DMRG calculations were found to be very close to the 
corresponding FCI results for the two small molecules, 
we found that the effect of truncation of the active 
space is smaller for DMRG-srPBE than for standard 


DMRG calculations. This effect was also visible in the 
entanglement-entropy measures considered. We studied 
the effect of the simultaneous treatment of static and 
dynamic correlation on such orbital-based descriptors, 
namely on single-orbital entropies and mutual informa¬ 
tion. We find that (i) the single-orbital entropies and 
mutual information are consistently smaller for DMRG- 
srDFT than for DMRG and that (ii) the major part of 
dynamical correlation is assigned to the short-range DFT 
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part so that the pattern of these entropy measures hardly 
change when the active space is reduced in a DMRG- 
srDFT calculation. 

The discussion of ligand-binding energies revealed the 
true potential of the DMRG-srDFT approach. For two 
reactions out of the WCCR10 benchmark set of ligand¬ 
binding energies, we found that DMRG-srDFT yields a 
much less pronounced dependence of the reaction energy 
on the size of the active space such that more consistent 
results were obtained. Focusing on one functional (PBE 
and its short-range variant), we found for reaction 1 a 
significant improvement on both pure DMRG and pure 
PBE results. However, the situation was more delicate 
in the case of reaction 2 for which we found a dramatic 
dependence of the reaction energy on the inclusion of 
specific orbitals in the active space such that the pure 
DMRG result deviates significantly from the experimen¬ 
tal reference, indicating a very pronounced contribution 
of orbitals beyond the chosen active spaces to the corre¬ 
lation energy. Still, the final DMRG-srPBE result turned 
out to be very similar to the pure PBE result. However, 
we also noted that all data still deviated much from the 
experimental reference, which might even be taken as an 
indication to reinvestigate the accuracy of the experimen¬ 
tal value. 

With DMRG as the wave function part, DMRG- 
srDFT provides access to much larger complete active 
orbital spaces than those feasible with any traditional 
CAS-type approach combined with DFT. Hence, includ¬ 
ing dynamical correlation through a short-range density 
functional is a viable option to preserve the efficiency of 
DMRG calculations by avoiding standard perturbation- 
theory-based approaches. This facilitates calculations 
with long-range CAS-type wave functions such that all 
remaining approximations are buried in the approximate 
short-range DFT functional. 

Although we obtained encouraging results for our case 
studies when comparing our new DMRG-srDFT ap¬ 
proach with truncated active orbital spaces to stan¬ 
dard DMRG and FCI, further improvement is possible 
and work along these lines is in progress in our lab¬ 
oratory. Apart from an extension of our implementa¬ 
tion to a spin-unrestricted framework in the srDFT part, 
also spin-state-specific short-range functionals will fur¬ 
ther be crucial for molecules with an open-shell electronic 
structural. 

Alternative approaches for a simultaneous treatment 
of static and dynamic correlation in a hybrid DMRG ap¬ 
proach that avoid a range-separation ansatz exist. In 
future work, we are considering the pair-density func¬ 
tional theory which was recently put forward for MC- 
SCF method^! as we ll as the site occupation density 
functional theory proposed by Fromager 511 -. 

It should finally be emphasized that hybrids between 
DFT and wave function methods are also expected to 
have less dramatic dependence on the one-electron ba¬ 
sis set than standard wave function methods. We are 
currently investigating this more quantitatively. 


ACKNOWLEDGMENTS 

EDH thanks the Villum foundation for a post-doctoral 
fellowship. This work has been financially supported 
by ETH Zurich and the Schweizer Nationalfonds (SNF 
project 200021L_156598). HJAaJ thanks the Danish 
Council for Independent Research | Natural Sciences for 
financial support (grant number 12-127074). 

1 J. Olsen, Int. J. Quantum Chem. Ill, 3267 (2011). 

2 S. R. White, Phys. Rev. Lett. 69, 2863 (1992). 

3 S. R. White, Phys. Rev. B 48, 10345 (1993). 

4 U. SCHOLLWOCK, Ann. Phys. 326, 96 (2011). 

5 D. Yaron, E. E. Moore, Z. Shuai, and J. L. Bredas, J. Chem. 
Phys. 108, 7451 (1998). 

6 Z. Shuai, J. L. Bredas, A. Saxena, and A. R. Bishop, J. Chem. 
Phys. 109, 2549 (1998). 

7 G. Fano, F. Ortolani, and L. ZlOSl, J. Chem. Phys. 108, 9246 
(1998). 

8 S. R. White and R. L. Martin, J. Chem. Phys. 110, 4127 
(1999). 

9 S. Daul, I. Ciofini, C. Daul, and S. R. White, Int. J. Quantum 
Chem. 79, 331 (2000). 

10 A. O. Mitrushenkov, G. Fano, F. Ortolani, R. Linguerri, 
and P. Palmieri, J. Chem. Phys. 115, 6815 (2001). 
n G. K. L. Chan and M. Head-Gordon, J. Chem. Phys. 116, 
4462 (2002). 

12 6. Legeza and J. Solyom, Phys. Rev. B 68, 195116 (2003). 
13 G. Moritz and M. Reiher, J. Chem. Phys. 126, 244109 (2007). 
14 J. Rissler, R. M. Noack, and S. R. White, Chem. Phys. 323, 
519 (2006). 

15 A. O. Mitrushenkov, G. Fano, R. Linguerri, and P. Palmieri, 
Int. J. Quantum Chem. 112, 1606 (2012). 

16 S. Sharma, G. K. L. Chan, and M. Head-Gordon, J. Chem. 
Phys. 136, 124121 (2012). 

17 N. Nakatani and G. K.-L. Chan, J. Chem. Phys. 138, 134113 
(2013). 

18 D. Zgid and M. Nooijen, J. Chem. Phys. 128, 144115 (2008). 
19 D. Zgid and M. Nooijen, J. Chem. Phys. 128, 014107 (2008). 
20 Y. Kurashige and T. Yanai, J. Chem. Phys. 130, 234114 
(2009). 

21 Y. Kurashige, G. K. L. Chan, and T. Yanai, Nature Chem. 5, 
660 (2013). 

22 S. Wouters, P. A. Limacher, D. van Neck, and P. W. Ayers, 
J. Chem. Phys. 136, 134110 (2012). 

23 S. Wouters, W. Poelmans, P. W. Ayers, and D. van Neck, 
Comput. Phys. Commun. 185, 1501 (2014). 

24 S. Knecht, O. Legeza, and M. Reiher, J. Chem. Phys. 140, 
041101 (2014). 

25 S. F. Keller and M. Reiher, Chimia 68, 200 (2014). 

26 S. Keller, M. Dolfi, M. Reiher, and M. M. Troyer, QC- 
MAQUIS: A massively parallel, Matrix-Product-Operator based 
Quantum-Chemical DMRG program., in preparation, 2015. 

27 T. Dresselhaus, J. Neugebauer, S. Knecht, S. F. Keller, 
Y. Ma, and M. Reiher, J. Chem. Phys 142, 044111 (2015). 
28 0. Legeza, R. Noack, J. Solyom, and L. Tincani, Applica¬ 
tions of Quantum Information in the Density-Matrix Renormal¬ 
ization Group, in Computational Many-Particle Physics, edited 
by H. Fehske, R. Schneider, and A. Weisse, volume 739 of 
Led. Notes Phys., pp. 653-664, Springer, Berlin/Heidelerg, 2008. 
29 G. K. L. Chan, J. J. Dorando, D. Ghosh, J. Hachmann, 
E. Neuscamman, H. Wang, and T. Yanai, Frontiers in Quan¬ 
tum Systems in Chemistry and Physics, chapter “An Introduc¬ 
tion to the Density Matrix Renormalization Group Ansatz in 
Quantum Chemistry”, p. 49, Springer, 2008. 

30 K. H. Marti and M. Reiher, Z. Phys. Chem. 224, 583 (2010). 
31 G. K. L. Chan and S. Sharma, Annu. Rev. Phys. Chem. 62, 
465 (2011). 


13 


32 K. H. Marti and M. Reiher, Phys. Chem. Chem. Phys. 13, 
6750 (2011). 

33 S. Wouters and D. van Neck, Eur. Phys. J. D 68, 272 (2014). 
34 Y. Kurashige, Mol. Phys. 112, 1485 (2014). 

35 S. Szalay, M. Pfeffer, V. Murg, G. Barcza, F. Verstraete, 
and O. Legeza, Tensor product methods and entanglement op¬ 
timization for ab initio quantum chemistry, arXiv:1412.5829vl, 
2014. 

36 M. Saitow, Y. Kurashige, and T. Yanai, J. Chem. Phys. 139, 
044118 (2013). 

37 Y. Kurashige, J. Chalupsky, T. Nguyen Lan, and T. Yanai, 
J. Chem. Phys. 141, 174111 (2014). 

38 T. Yanai, Y. Kurashige, W. Mizukami, J. Chalupsky, 
T. Nguyen Lan, and M. Saitow, Int. J. Quantum Chem. 115, 
283 (2015). 

39 I. Shavitt, Int. J. Mol. Sci. 3, 639 (2002). 

40 S. Sharma and G. K.-L. Chan, J. Chem. Phys. 141, (2014). 
41 J. P. Perdew, A. Savin, and K. Burke, Phys. Rev. A 51, 4531 
(1995). 

42 J. M. Seminario, Int. J. Quantum Chem. S28, 655 (1994). 

43 J. P. Perdew and K. Burke, Int. J. Quantum Chem. 57, 309 
(1996). 

44 0. Gunnarson and B. I. Lundqvist, Phys. Rev. B 13, 4274 
(1976). 

45 A. Savin and H.-J. Flad, Int. J. Quantum Chem. 56, 327 
(1995). 

46 S. Grimme and M. Waletzke, J. Chem. Phys. Ill, 5645 (1999). 
47 A. Savin, Recent Developments and Applications of Modern 
Density Functional Theory , p. 327, Elsevier, Amsterdam, 1996. 
48 E. Fromager, J. Toulouse, and H. J. Aa. Jensen, J. Chem. 
Phys. 126, 074111 (2007). 

49 E. Fromager, R. Cimiraglia, and H. .J. Aa. Jensen, Phys. 
Rev. A 81, 024502 (2010). 

50 G. L. Manni, R. K. Carlson, S. Luo, D. Ma, J. Olsen, D. G. 
Truhlar, and L. Gagliardi, J. Chem. Theory and Comput. 
10, 3669 (2014). 

51 E. Fromager, Mol. Phys. 113, 419 (2015). 

52 C. M. Marian and N. Gilka, J. Chem. Theory Comput. 4. 1501 
(2008). 

53 M. R. Silva-Junior, M. Sohreiber, S. P. A. Sauer, and 
W. Thiel, J. Comput. Phys. 129, 104103 (2008). 

54 M. M. Alam and E. Fromager, Chem. Phys. Lett. 554, 37 

( 2012 ). 

55 E. Fromager, S. Knecht, and H. .J. Aa. .Jensen, J. Chem. 
Phys. 138, 084101 (2013). 

56 E. D. Hedegard, N. H. List, H. J. Aa. Jensen, and J. Kong- 
STED, J. Chem. Phys. 139, 044101 (2013). 

57 V. Veryazov, P.-A. Malmqvist, and B. O. Roos, Int. J. Quan¬ 
tum Chem. Ill, 3329 (2011). 

58 L. M. Cheung, K. R. Sundberg, and K. Ruedenberg, Int. J. 

Quantum Chem. 16. 1103 (1979). 

59 H. J. Aa. Jensen, P. Jorgensen, H. Agren, and J. Olsen, J. 
Chem. Phys. 88, 3834 (1988). 

60 P. Pulay and T. P. Hamilton, J. Chem. Phys. 88, 4926 (1988). 
61 6. Legeza and J. Solyom, Phys. Rev. Lett. 116401, 116401 
(2006). 

62 K. BoguSLAWSKI and P. Tecmer, Int. J. Quantum Chem. 
(2014). 

63 K. Boguslawski, P. Tecmer, 6. Legeza, and M. Reiher, J. 

Phys. Chem. Lett. 3, 3129 (2012). 

64 T. Weymuth, E. P. A. Couzijn, P. Chen, and M. Reiher, J. 

Chem. Theory and Comput. 10, 3092 (2014). 

65 T. Weymuth and M. Reiher, Int. J. Quantum Chem. 115, 90 
(2015). 

66 B. O. ROOS, P. R. Taylor, and P. E. M. Siegbahn, Chem. 
Phys. 48, 157 (1980). 

67 P. E. M. Siegbahn, J. Almlof, A. Heiberg, and B. O. Roos, 
J. Chem. Phys. 74, 2384 (1981). 

68 T. Helgaker, P. Jorgensen, and J. Olsen, Molecular 
Electronic-Structure Theory , Wiley (pages 482 and 645-646), 


2004. 

69 E. Fromager and H. J. Aa. Jensen, Phys. Rev. A 78, 022504 
(2008). 

70 E. D. Hedegard, F. Heiden, S. Knecht, E. Fromager, and 
H. J. Aa. Jensen, J. Chem. Phys. 139, 184308 (2013). 

71 E. Goll, H.-J. Werner, and H. Stoll, Phys. Chem. Chem. 
Phys. 7, 3917 (2005). 

72 P. Katarzyna, J. Chem. Phys. 136, 184105 (2012). 

73 H. Ikura, T. Tsuneda, T. Yanai, and K. Hirao, J. Chem. 
Phys. 115, 3540 (2001). 

74 Y. Tawada, T. Tsuneda, S. Yanagisawa, T. Yanai, and K. Hi¬ 
rao, J. Chem. Phys. 120, 8425 (2004). 

75 0. A. Vydrov, J. Heyd, A. V. Krukau, and G. E. Scuseria, 

J. Chem. Phys. 125, 074106 (2006). 

7e O. A. Vydrov and G. E. Scuseria, J. Chem. Phys. 125, 234109 
(2006). 

77 I. C. Gerber and J. G. Angyan, Chem. Phys. Lett. 415, 100 
(2005). 

78 J. K. Pedersen, Description of Correlation and Relativistic Ef¬ 
fects in Calculations of Molecular Properties, PhD thesis, Uni¬ 
versity of Southern Denmark, 2004. 

79 Z. Huang and S. Kais, Chem. Phys. Lett. 413, 1 (2005). 
80 DALTON, a molecular electronic-structure program, Release 
Dalton2015 (2015), see http://daltonprogram.org/ development 
version, 2015. 

81 K. Aidas, C. Angeli, K. L. Bak, V. Bakken, R. Bast, L. Bo- 
man, O. Christiansen, R. Cimiraglia, S. Coriani, P. Dahle, 
E. K. Dalskov, U. Ekstrm, T. Enevoldsen, J. J. Erik- 
sen, P. Ettenhuber, B. Fernandez, L. Ferrighi, H. Fliegl, 
L. Frediani, K. Hald, A. Halkier, C. Hattig, H. Heiberg, 
T. Helgaker, A. C. Hennum, H. Hettema, E. HjertenjES, 

S. Host, I.-M. Hoyvik, M. F. Iozzi, B. Jansi'k, H. J. Aa. 
Jensen, D. Jonsson, P. Jorgensen, J. Kauczor, S. Kirpekar, 

T. Kjjergaard, W. Klopper, S. Knecht, R. Kobayashi, 
H. Koch, J. Kongsted, A. Krapp, K. Kristensen, A. Ligabue, 
O. B. LUTNiES, J. I. Melo, K. V. Mikkelsen, R. H. Myhre, 
C. Neiss, C. B. Nielsen, P. Norman, J. Olsen, .J. M. H. 
Olsen, A. Osted, M. J. Packer, F. Pawlowski, T. B. Ped¬ 
ersen, P. F. Provasi, S. Reine, Z. Rinkevicius, T. A. Ruden, 

K. Ruud, V. V. Rybkin, P. Salek, C. C. M. Samson, A. S. 
de Meras, T. Saue, S. P. A. Sauer, B. Schimmelpfennig, 
K. Sneskov, A. H. Steindal, K. O. Sylvester-Hvid, P. R. 
Taylor, A. M. Teale, E. I. Tellgren, D. P. Tew, A. J. Thor¬ 
valdsen, L. Thogersen, O. Vahtras, M. A. Watson, D. J. D. 
Wilson, M. Ziolkowski, and H. Agren, WIREs Comput. Mol. 
Sci. 4, 269 (2013). 

82 T. H. Dunning Jr., J. Chem. Phys. 90, 1007 (1989). 

83 J. Heyd, G. E. Scuseria, and M. M. Ernzerhof, J. Chem. 
Phys. 118, 8207 (2015). 

84 J. Olsen, P. Jorgensen, H. Koch, A. Balkova, and R. J. 

Bartlett, J. Chem. Phys. 104, 8007 (1996). 

85 G. K. L. Chan, M. Kallay, and J. Gauss, J. Chem. Phys. 121, 
6110 (2004). 

86 A. D. Becke, Phys. Rev. A 38, 3098 (1988). 

87 F. Weigend and R. Ahlrichs, Phys. Chem. Chem. Phys. 7, 
3297 (2005). 

88 F. Weigend and R. Ahlrichs, Phys. Chem. Chem. Phys. 8, 
1057 (2006). 

89 R. Ahlrichs, M. Bar, M. Haser, H. Horn, and C. Kolmel, 
Chem. Phys. Lett. 162, 165 (1989). 

90 J. P. Perdew, K. Burke, and M. Ernzerhof, Phys. Rev. Lett. 
77, 3865 (1996). 

91 A. Schafer, C. Huber, and R. Ahlrichs, J. Chem. Phys. 100, 
5829 (1994). 

92 See supplemental material at [URL will be inserted by AIP] for 
a discussion of total electronic energies for H 2 O and N 2 and for 
graphical representations of orbitals of the active space in DMRG 
calculations on the ligand-dissociation reactions. 

93 C. R. Jacob and M. Reiher, Int. J. Quantum Chem. 112, 3661 

( 2012 ). 


Electronic Supporting Infomation: 


Density Matrix Renormalization Group with Efficient Dynamical Electron Correlation 

Through Range Separation 

Erik Donovan Hedegard, 1, a ) Stefan Knecht, 1 Jesper Skau Kielberg, 2 Hans J0rgen 
Aagaard Jensen, 2,b ) and Markus Reiher 1,c ) 

^ETH Zurich, Laboratorium fur Physikalische Chemie, Vladimir-Prelog Weg-2, 
CH-8093 Zurich, Switzerland 

2) Department of Physics, Chemistry and Pharmacy, University of Southern Denmark, 
Campusvej 55, Odense, Denmark 

(Dated: May 27, 2015) 


Electronic mail: erik.hedegard@phys.chem.ethz.ch 
b 'Electronic mail: hjj@sdu.dk 

c )Electronic mail: markus.reiher@phys.chem.ethz.ch 


1 



I. TOTAL ELECTRONIC ENERGIES FOR WATER AND DINITROGEN 


H 2 0 and N 2 have been investigated extensively with both traditional wave function meth¬ 
ods and DMRG 1 3 . Here, we compare to the work by Olsen et al 4 for H 2 0 and the studies 
by Krogh and Olsen 5 , Larsen et al. 6 , and Chan et al. 2 for N 2 . 

In the following, the DMRG and DMRG-srPBE results are compared in terms of the to¬ 
tal energy and of the non-parallelity error (NPE). The NPE measures the difference between 
the absolute maximum and minimum deviations from the DMRG-FCI potential curve over 
the internuclear distances considered (see e.g. Ref. 7). The NPE would be zero if the curves 
differed by a constant shift. As the NPE suffers from that the individual methods do not 
share the same energy scale (due to the approximate srDFT functional) we have addition¬ 
ally reported the energy difference (in kJ mol -1 ) with respect to the electronic equilibrium 
internuclear distance (R e ). 

Starting with H 2 0, the energies calculated at several distances along the symmetric 
stretching coordinate are summarized in Table I. The DMRG(10,24)[512] data agree well 
with the FCI results by Olsen et al. 4 which also shows that the number of renormalized active- 
subsystem block states m — 512 is sufficient to achieve FCI accuracy with DMRG. It is fur¬ 
ther seen that the DMRG-srPBE energies are higher than the energies of DMRG( 10,24) [512], 
even for DMRG(10,24)[512]-srPBE where a FCI long-range wave function is employed. Al¬ 
though one might intuitively expect that additional dynamical correlation would lower the 
energy, the approximate nature of the srDFT functionals should again not be forgotten. As 
mentioned in the Theory section of the main article, there is no guarantee that the WFT- 
srDFT energy for a given ^-value is below the corresponding WFT energy. A similar result 
has been observed by Fromager et al. 8 

Finally, we should note that the DMRG(10,24)[512]-srPBE data can be considered as 
a technical consistency test of the approach. A CAS(10,24) active space represents the 
full configuration interaction scenario, for which the srDFT should not give any additional 
contribution within the one-electron basis set chosen. For H 2 0 we understand from the 
DMRG(10,24)[512]-srPBE result in Table I that 918 — 762 = 156 kJ/mol are artificially 
introduced due to two technical limitations: (i) the spin-restricted implementation of the 
srDFT part and (ii) the approximate nature of the short-range functional. 

Turning to N 2 , all results are compiled in Table II. Our reference DMRG(14,28)[2048] 
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results are in good agreement with the results reported by Chan and co-workers 2 , although 
small differences occur far from the equilibrium distance, presumably due to different pro¬ 
cedures during the DMRG start-up sweeps. 

For both H 2 0 and N 2 , the effect of truncation of the active space is smaller for DMRG- 
srPBE than for standard DMRG. In the main paper, we have elaborated on this issue by 
analysis of the entanglement entropies. 

With respect to the NPEs of DMRG-srPBE, Tables I and II take the standard FCI data 
as reference. For this reference, the NPEs are for H 2 0 0.0595 (DMRG(10,24)[512]-srPBE) 
and 0.0592 (DMRG(8,8)[512]-srPBE), respectively. Thus, for H 2 0, the NPE of DMRG- 
srPBE is in between the NPEs of PBE (0.0342) and PBE0 (0.0677). The NPEs are for N 2 
0.0986 (DMRG(14,28)[1024]-srPBE) and 0.0984 (DMRG(6,6)[1024]-srPBE), whereas the 
PBE functional gives an NPE of 0.0691. We note that for N 2 , the PBE0 functional even 
crosses the DMRG(14,28)[2048] curve, which is neither observed for DMRG-srPBE nor 
for standard PBE, and which makes it impossible to calculate a well-defined NPE. More 
importantly, the NPEs of the truncated-active-space DMRG-srPBE are in both cases, H 2 0 
and N 2 , nearly identical up to A(NPE) = 0.0003 with their DMRG(FCI)-srPBE references. 
This substantiates the above finding that even for truncated active spaces accurate results 
close to DMRG(FCI)-srDFT quality are within reach; an observation that will be important 
for larger molecules where FCI is beyond the limits, even for DMRG. 

It should also be mentioned that the NPEs for both H 2 0 and N 2 with DMRG-srDFT 
are larger than for pure DMRG, also when using a DMRG(FCI)-srDFT description. The 
same holds for the energies relative to R e , which are overestimated (this is also the case for 
regular DFT). Since one for DMRG(FCI)-srDFT have an exact description of the long-range 
wave function, the errors arise from the approximate srDFT functional. One obvious source 
of error is that the used functionals are not spin-adapted. Further, the applied short-range 
functional is not self-interaction free, as discussed in Ref. 8. 

From the discussion of the total electronic energies we can conclude that: 

• The small active space is already so large for these small molecules that dynamical 
correlation plays a minor role. Consequently, relative energies are already very well 
reproduced by the small CAS. 

• Still, it can be seen that the effect of truncation of the active space is smaller for 
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TABLE I. DMRG, DMRG-srDFT, FCI, and DFT energies for H 2 0 at various O—H distances 
along the symmetric streching mode in a Dunning cc-pVDZ basis set. Values in parantheses are 
differences (in kJ moL 1 ) with respect to the electronic equilibrium internuclear distance (i? e ) for a 
given setup. DMRG-srDFT data at increasing distance become unreliable as our implementation 
is a spin-restructed one that cannot properly describe the open-shell fragments, ao denotes the 
Bohr radius. 


Method 1.843 ao (Re) 2.765 ao 3.687 a 0 NPE 


FCI = CAS(10,24) (Ref. 4) 

-76.241 

860 

(0.0) - 

-76.072 

348 

(445) - 

-75.951 

711 

(762) 

— 

DMRG(10,24) [512] 

-76.241 

748 

(0.0) - 

-76.072 

206 

(445) - 

-75.951 

577 

(762) 

— 

DMRG(8,8)[512] 

-76.064 

887 

(0.0) - 

-75.906 

261 

(417) - 

-75.797 

043 

(703) 

0.0223 

DMRG(10,24) [512]-srPBE 

-76.334 

180 

(0.0) - 

-76.153 

803 

(474) - 

-75.984 

516 

(918) 

0.0595 

DMRG(8,8) [512]-srPBE 

-76.331 

452 

(0.0) - 

-76.150 

984 

(474) - 

-75.982 

104 

(917) 

0.0592 

PBE 

-76.333 

974 

(0.0) - 

-76.173 

674 

(421) - 

-76.018 

828 

(827) 

0.0342 

PBEO 

-76.338 

765 

(0.0) - 

-76.155 

765 

(481) - 

-75.980 

886 

(940) 

0.0677 


DMRG-srPBE than for standard DMRG. 

• Total electronic energies at increasing internuclear distance and thus relative energies 
are unreliable because of a technical limitation of our current implementation (that is a 
restriction to spin-restricted electronic structures, which is not valid in the asymptotic 
limit). 


II. LIGAND-BINDING ENERGIES 


To demonstrate that convergence w.r.t. the DMRG macro-iterations (sweeps) has been 
achieved, we present the corresponding data in Tables III and IV. Total electronic energies 
for the ligand-dissociation reactions are reported in Tables V and VI. 
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TABLE III. Dissociation energies for reaction 1, obtained with different number of renormalized 
states (m) and sweeps. The results have not been corrected for the zero-point vibrational energy. 


Method 

Sweeps 

D e (kJ/mol) 

DMRG[2000] (30,22) 

20 

173.5 

DMRG[2000] (30,22) 

10 

173.5 

DMRG[1000] (30,22) 

10 

173.5 

DMRG[2000] (30,22)-srPBE 

20 

225.1 

DMRG[2000] (30,22)-srPBE 

10 

225.1 

DMRG[1000](30,22)-srPBE 

10 

225.1 


TABLE IV. Dissociation energies for reaction 2, obtained with different number of renormalized 
states (m) and sweeps. The results have not been corrected for the zero-point vibrational energy. 


Method 

Sweeps 

D e (kJ/mol) 

DMRG[2000] (22,20) 

20 

34.0 

DMRG[2000] (22,20) 

10 

34.0 

DMRG[1000] (22,20) 

10 

34.0 

DMRG[2000] (22,20)-srPBE 

20 

81.6 

DMRG[2000] (22,20)-srPBE 

10 

81.6 

DMRG[1000](22,20)-srPBE 

10 

81.7 


TABLE V. Total electronic energies in Hartree for ligand-dissociation reaction 1 along with corre¬ 
sponding dissociation energies (def-TZVP basis set). Dq is zero-point vibrational-energy corrected 
(see the main paper). 


Method 

Reactant 

Product D e (kJ/mol) D 0 (kJ/mol) 

DMRG[2000] (30,22) 

-2129.899 

308 

-2129.832 

571 

173.5 

165.1 

DMRG[2000] (20,18) 

-2129.896 

908 

-2129.831 

567 

169.9 

161.5 

DMRG[2000] (10,10) 

-2129.880 

839 

-2129.829 

791 

132.8 

124.3 

DMRG[2000] (30,22)-srPBE 

-2133.524 

563 

-2133.438 

019 

225.1 

216.6 

DMRG[2000](20,18)-srPBE 

-2133.524 

455 

-2133.436 

822 

227.9 

219.4 

DMRG[2000](10,10)-srPBE 

-2133.518 

977 

-2133.435 

733 

216.5 

208.0 

DMRG[2000](30,22)-srPBE(GWS) -2133.609 491 

-2133.514 696 

246.5 

238.1 
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TABLE VI. Total electronic energies in Hartree for ligand-dissociation reaction 2 along with corre¬ 
sponding dissociation energies (def-TZVP basis set). Dq is zero-point vibrational-energy corrected 
(see the main paper). 


Method 

Reactant 

Product 

D e (kJ/mol) 

D 0 (kJ/mol) 

DMRG [2000] (24,24) 

-712.280 

859 

-712.263 842 

44.3 

36.9 

DMRG[2000] (22,20) 

-712.269 

411 

-712.256 329 

34.0 

26.7 

DMRG[2000] (18,18) 

-712.265 

749 

-712.251 309 

37.6 

30.2 

DMRG[2000] (16,16) 

-712.260 

395 

-712.238 611 

56.6 

49.3 

DMRG[2000] (14,14) 

-712.255 

088 

-712.233 875 

55.2 

47.8 

DMRG[2000] (10,10) 

-712.251 

167 

-712.227 939 

60.4 

53.1 

DMRG [2000] (8,8) 

-712.249 

671 

-712.224 544 

65.3 

58.0 

DMRG [2000] (22,20)-srPBE 

-716.065 

249 

-716.0338 492 

81.7 

74.3 

DMRG[2000](8,8)-srPBE 

-716.058 

103 

-716.0271 214 

80.6 

73.2 

DMRG [2000] (22,20)-srPBE(GWS) 

-716.065 

249 

-716.033 849 

101.8 

94.5 


TABLE VII. Natural orbital occupation numbers for reactant and product structures of reaction 


2. 

The orbitals are 

shown in 

Figures 5-8. 






Orbital Reactant 

Product 

Reactant 

Product 

Reactant 

Product 

Reactant 

Product 


DMRG(16,16) 

DMRG(18,18) 

DMRG(22,20) 

DMRG(22,20)-srPBE 

54 

Inactive 

Inactive 

Inactive 

Inactive 

1.99927 

1.99999. 

1.99975 

1.99984 

55 

Inactive 

Inactive 

Inactive 

Inactive 

1.99983 

1.98724 

1.99982 

1.99053 

56 

Inactive 

Inactive 

1.99965 

1.96652 

1.99947 

1.96575 

1.99928 

1.99420 

57 

1.99248 

1.99860 

1.99393 

1.99857 

1.99521 

1.99851 

1.99986 

1.99999. 

58 

1.99927 

1.99892 

1.99785 

1.99950 

1.99665 

1.99950 

1.99214 

1.99621 

59 

1.99979 

1.99618 

1.99958 

1.99571 

1.99959 

1.99868 

1.99658 

1.99774 

60 

1.99874 

1.97459 

1.99822 

1.97403 

1.99303 

1.96981 

1.99984 

1.99642 

61 

1.99515 

1.99124 

1.99166 

1.99003 

1.98939 

1.98987 

1.99423 

1.99999. 

62 

1.97704 

1.99999. 

1.97567 

1.99999. 

1.97617 

1.99999. 

1.99581 

1.99999. 

63 

1.99611 

1.99999. 

1.99002 

1.99999. 

1.99056 

1.99999. 

1.99914 

1.99916 

64 

1.99787 

1.99224 

1.99615 

1.99265 

1.99614 

1.99273 

1.99560 

1.99987 

65 

0.02884 

0.02590 

0.03046 

0.02689 

0.03189 

0.03606 

0.01009 

0.00597 

66 

0.00005 

0.00638 

0.00005 

0.03751 

0.00008 

0.03861 

0.00675 

0.01128 

67 

0.00002 

0.00001 

0.00003 

0.00001 

0.00006 

0.00003 

0.00631 

0.00646 

68 

0.00493 

0.00005 

0.01314 

0.00003 

0.01347 

0.00004 

0.00003 

0.00001 

69 

0.00761 

0.00047 

0.00902 

0.00071 

0.00983 

0.00075 

0.00066 

0.00175 

70 

0.00032 

0.01288 

0.00010 

0.01470 

0.00013 

0.01512 

0.00005 

0.00045 

71 

0.00012 

0.00004 

0.00017 

0.00010 

0.00023 

0.00013 

0.00065 

0.00005 

72 

0.00165 

0.00251 

0.00179 

0.00295 

0.00467 

0.00705 

0.00333 

0.00006 

73 

Secondary 

Secondary 

0.00250 

0.00011 

0.00432 

0.00012 

0.00006 

0.00003 

74 

Secondary 

Secondary 

Secondary 

Secondary 

Secondary 

Secondary 

Secondary 

Secondary 
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DMRG[2000](30,22) Active Space Orbitals 
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FIG. 1. Active space orbitals for the complex in reaction 1. The occupied orbitals are underlined 


in bold font. 
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DMRG[2000](30,22) Active Space Orbitals 
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FIG. 2. Active space orbitals for the stretched structure (Cu-N internuclear distance 7 A) in 
reaction 1. The occupied orbitals are underlined in bold font. 
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DMRG[2000](30,22)-srPBE Active Space Orbitals 


^ A* 

42 43 44 45 

I** %. & %. 

46 47 48 49 



62 63 


FIG. 3. Active space orbitals for the complex in reaction 1. The occupied orbitals are underlined 
in bold font. 
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DMRG[2000](30,22)-srPBE Active Space Orbitals 
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FIG. 4. Active space orbitals for the stretched structure (Cu-N internuclear distance 7 A) in 
reaction 1. The occupied orbitals are underlined in bold font. 
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DMRG[2000](22,20) Active Space Orbitals 
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FIG. 5. Active space orbitals for the complex in reaction 2. The occupied orbitals are underlined 
in bold font. 
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DMRG[2000](22,20) Active Space Orbitals 
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FIG. 6. Active space orbitals for the stretched structure (Pt-N internuclear distance 7 A) in 
reaction 2. The occupied orbitals are underlined in bold font. 
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DMRG[2000](22,20)-srPBE Active Space Orbitals 



FIG. 7. Active space orbitals for the complex in reaction 2. The occupied orbitals are underlined 
in bold font. 
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DMRG[2000](22,20)-srPBE Active Space Orbitals 



FIG. 8. Active space orbitals for the stretched structure (Pt-N internuclear distance 7 A) in 
reaction 2. The occupied orbitals are underlined in bold font. 
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